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Capillary forces significantly affect the stability of sandpiles. We analyze the stability of sandpiles 
with such forces, and find that the critical angle is unchanged in the limit of an infinitely large 
system; however, this angle is increased for finite-sized systems. The failure occurs in the bulk of 
the sandpile rather than at the surface. This is related to a standard result in soil mechanics. The 
increase in the critical angle is determined by the surface roughness of the particles, and exhibits 
three regimes as a function of the added-fluid volume. Our theory is in qualitative agreement with 
the recent experimental results of Hornbaker et al., although not with the interpretation they make 
of these results. 
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The continuum mechanics of most materials was es- 
tablished in the 19th century; however, the mechanics 
of granular materials is still largely mysterious [jjj . The 
study of granular media is also motivated by the ubiquity 
of this form of matter in a variety of industrial contexts, 
as well as in geophysical ones. 

While most recent attention has focussed on dry gran- 
ular media, a recent experimental study by Hornbaker et 
al. has opened the relatively unexplored subject of "hu- 
mid" granular media, in which small amounts of added 
fluid generate, through capillarity, adhesive forces be- 
tween the grains Somewhat whimsically, these au- 
thors argue that their work is appropriate for the un- 
derstanding of sandcastles; we actually take this point 
seriously, because adhesive forces and other liquid effects 
are extremely important in geophysical applications, of 
which sandcastles are an unusual example. 

In this Letter, we present a theory of the stability of hu- 
mid sandpiles, based upon a continuum analysis of their 
statics. While previous work on the statics and dynam- 
ics of dry sandpiles concentrated on the behavior of the 
pile's surface ||, we find that the addition of small ad- 
hesive forces between the grains causes the site of failure 
to move from the surface into the bulk of the sandpile, a 
fact well-known in soil mechanics. Even though the fail- 
ure of the sandpile at the critical angle occurs at some 
finite depth, in the limit of infinite system size, the crit- 
ical angle is actually unchanged by the adhesion. For 
finite systems the angle of repose is increased from the 
infinite-system/non- adhesive critical angle. By analyzing 
the cohesive effect of small amounts of wetting fluid, we 
find that this increase in the critical angle as a function 
of the added-fluid volume exhibits a range over which the 
dependence is linear, in agreement with the principal re- 
sult of Hornbaker et al. However, we disagree with the 
suggestion of these authors that most of the wetting fluid 
will be found outside the particle contact zones. 

To determine the stability of a sandpile, we must have 
a criterion for local failure of the sandpile. For a non- 
adhesive (dry) sandpile, a simple phenomenological cri- 
terion for failure is that 



r > ka, (1) 

where r is the tangential stress across some plane interior 
to the sandpile, a is the normal compressive stress across 
that plane, and k is the internal friction coefficient. For a 
given stress state stability requires that there be no plane 
for which the ratio of r to a exceeds k. To determine the 
plane on which this ratio is maximized we turn to the 
Mohr construction ||. 

The Mohr circle provides a geometric construction to 
transform under rotations a two dimensional symmetric 
tensor such as the stress tensor. Suppose the two prin- 
cipal stresses at a point are <j\ and &i. Then a circle 
drawn in the a — t plane through the points (ci , 0) and 
(<72,0), symmetrically about the a axis, gives the nor- 
mal and tangential stresses on any plane. To find these 
stresses for a physical system rotated by an angle <fi from 
the principal system one looks at the points on the circle 
at an angle 2(f> from the horizontal. These two opposite 
points give the normal and tangential stresses across the 
physical planes of the rotated system. Thus the stability 
criterion for a dry sandpile is that no Mohr circle may 
extend above the line r = ka. 

Consider a semi-infinite dry sandpile, whose surface 
is oriented at an angle 9 to the horizontal. We choose 
an x — z coordinate system, in which z gives the distance 
from the surface of the pile (z > down) , and x gives the 
distance parallel to the surface. Then the stress tensor 
Oij satisfies the static equations 

d z cr zz + d x o- xz = pg cos 0, (2) 
d z a zx + d x cr xx = pg sin 9, (3) 

where p is the density of the sandpile. To solve these 
equations, we first restrict ourselves to solutions which 
are functions of z alone-in a semi-infinite system, any 
x-dependence of the solutions would be liable to gen- 
erate arbitrarily large stresses near the surface, which 
would cause the system to buckle. The most general 
z-dependent solution, which also satisfies the boundary 
condition that the surface is stress-free, is 
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g zz = pgz cos 0, 
<?xz = pgz sin 0, 
Oxx = C{z), 



(4) 
(5) 
(6) 



The well-known stress indeterminacy in granular media 
implies that C(z) is an unknown function. We can fix this 
function, however, by finding the form of C(z) that al- 
lows us to maximize the critical angle. The angle thereby 
obtained will certainly be an upper bound on the true 
critical angle; furthermore, if the sandpile is free to ad- 
just its undetermined stress C(z) within some range, we 
expect that this upper bound will be identical to the crit- 
ical angle. Thus, while it is impossible in general to fix 
C(z) without some constitutive or closure relation, we do 
believe that it is possible to fix it at the critical angle. 

Clearly the function C(z) that will maximize the crit- 
ical angle will have the form C(z) = C'z. To apply the 
geometric Mohr construction, we first determine the di- 
ameter and the position of the center of the Mohr cir- 
cle given by the difference and sum of eigenvalues of the 
stress tensor respectively. From Eqs. we find that 

these eigenvalues are 



01,2 = zpg cos 



1 



± 



(1 



tan 2 



(7) 



where c = C'/(pg cos 0). Note that both the radius and 
center position of the Mohr circle depend linearly on the 
depth z below the surface. The maximum value of t/ct 
as a function of c is 



[l + cf 



4(c - tan 2 9) 



(8) 



Now this quantity must be < k, so to find the maximum 
critical angle, we wish to minimize the right-hand side 
of Eq. (||) with respect to c. We then find that at this 
minimum value of c, 



tan 9, 



(9) 



so that the critical angle is 9 C = tan -1 k. which is a classic 
result 0]. This simple model, by itself, does not indicate 
at which depth failure initiates-presumably a dynamical 
model would resolve this ambiguity l|. 

Now consider a sandpile in which a normal adhesive 
stress sa is exerted across every plane, in addition to 
whatever other stresses may exist due to the body forces. 
This stress introduces a normal force between pairs of 
contiguous particles which allows the sandpile to support 
a finite shear stress, even in the limit of zero applied com- 
pressive stress. The maximum supported shear stress, in 
this case, is ksA and we therefore replace the dry sand 
failure criterion, Eq. ([!]), by 



t > k(a + sa) 



(10) 



Performing a calculation similar to the one above, we 
find that the failure criterion, which is now an explicit 
function of depth, is 



k = tan 0(1 



SA 



pgz cos 



(11) 



Note that since the sandpile will only be in a state of 
incipient failure below some fixed depth, the stress tensor 
retains some indeterminacy above that depth. Our Mohr 
analysis, now local, applies only at the incipent failure 
depth and does not determine the global stress state of 
the sand pile. 

The criterion thus derived, Eq. (|TT|) , is most stringent 
as z — > oo, in which case the dry sandpile result C = 
tan -1 k is recovered. On the other hand, for a sandpile 
of fixed depth D, the failure must occur at most at depth 
D. Thus the critical angle will be the solution of the 
equation 



k = tan C (£>)(! 



pgD cos 9 C (D) ' 



(12) 



giving an critical angle 9 C (D) that decreases monotoni- 
cally with D. Thus finite humid sandpiles have a depth- 
dependent critical angle, unlike dry sandpiles, which is 
a well-known result in soil mechanics 0|. In addition, 
we see that humid (i.e. cohesive) sandpiles fail at depth, 
whereas from statics alone we were unable to determine 
the failure depth of dry sandpiles. In the case of small 
adhesion stress sa/ pgD < 1, we can write 



tan 9 r ~ k 



ks A 

pgD 



sec [tan 1 (fc)] 



(13) 



If, as we are assuming, the adhesion arises from cap- 
illary forces, we must still connect the adhesive stress to 
the amount of fluid present. We suppose that the sand 
is composed of macroscopically spherical grains (radius 
R) whose surface roughness may be characterized as fol- 
lows: the spatial correlation of fluctuations in local sur- 
face height saturates at height In at a lateral distance d 
that is much smaller than the particle radius, d <C R. 
Since the particles are macroscopically spherical, we re- 
quire that Ir <c R (see Fig. |l|). 

We can characterize the surface roughness of two par- 
ticles in contact by considering the function 5(x) which 
gives the average distance between the two particles a 
lateral distance x from an asperity at which the two par- 
ticles are in contact. We write S(x) in the form Q 



5(x) = l R f(x/d) , 
where f(w) is a scaling function with the limits 



w 
w 



0. 

oo . 



(14) 



(15) 



The roughness exponent, x, of the surface satisfies < 
X<1- 
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FIG. 1. The contact zone between two rough particles of 
radius R. The scale of height deviations from the mean is Ir 
and the height fluctuations are correlated over a distance d. 
The lateral size of the contact zone in which the macroscopic 
curvature of the particles is not apparent is ~ \/IrR. 

Note that this form can only be valid for x < VTrR. 
For larger values of x, the macroscopic curvature of the 
particles will determine the local distance between them 
(see Fig. |). 

For such particles, there are three regimes for the cap- 
illary force exerted by a wetting fluid as a function of V, 
the total amount of fluid present per particle contact. 

Asperity Regime- For the smallest values of V, the 
capillary force is dominated by the accumulation of fluid 
around a single or a small number of asperities at which 
two neighboring particles are in contact. This will hold 
until the lateral extent of the fluid-filled region exceeds 
d, determining a maximum contact fluid volume for this 
regime, V\ ~ IrcP. 

We can write the adhesion force Ja as 



f - TA 

J A — , 

r 



(16) 



where T is the surface tension of the fluid, r is the radius 
of curvature of the meniscus of the fluid layer connect- 
ing the two grains near the asperity and A ~ V/5 is the 
area of the contact patch. —Y/r is the pressure reduc- 
tion due to the capillary meniscus. Because r will be 
approximately the distance between the particles at the 
meniscus, we find 



Ir \Vi 
and the adhesion force is 



(17) 
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for V < Vt, 



(18) 



where V\ = IrcP. 

For a rough surface where X = 1, Eq. (fFJJ) shows that 
the force depends on the cube root of the fluid volume. 
This is identical to the dependence of the cohesive force 
between a cone and plate on the volume of the liquid 
bridge connecting them || . It is to be expected that the 
cone-and-plate model will reproduce the cohesive force 
near a single asperity. 

Roughness Regime- For larger values of V, the fluid 
will occupy a statistically rough region, which is still 
small enough that the macroscopic curvature of the par- 
ticles plays no role-however, the fluid occupies more than 
the area around a single asperity. This regime occurs for 
Vi < V < V 2 , where V 2 



P 



l R R. The pressure is 

r 



and the force will be 



TV 

T 2 

l R 



f A ~ -g-, for Vi < V < V 2 



(19) 



(20) 



In this roughness regime, the cohesive force is linear in 
the volume of the added fluid, reproducing the linear 
dependence found by Hornbaker et al. 

Spherical Regime- When the lateral extent of the 
fluid contact exceeds d, then the wetting region will be 
determined by the macroscopic curvature of the particles, 
and the surface roughness will no longer play a significant 
role. In this case the pressure is given by M 



P = - 



r 



y/V/2wR ' 



and the force by 



f A = 2nTR, for V > Vi , 



(21) 



(22) 



which is independent of the volume of the liquid bridge 
joining the two grains. Thus the linear increase of the 
cohesive force with fluid volume saturates for volumes 
V > V 2 = l R R (see Fig. |). 

If the fluid wets the surface of the particles, then in 
addition to the fluid in the contact region, there will also 
be a layer of fluid of thickness t on the surface of the par- 
ticles. Typically this film will be no thicker than a few 
monolayers; hence there are a complicated set of forces 
between this film and the surface. To simplify, we con- 
sider only the Van der Waals forces, which for a thickness 
t generate a "disjoining pressure" Pa given by 



Pd = 



2H 



(23) 



where H is the Hamaker constant (H < for a wetting 
fluid). To determine the thickness of the wetting region, 



3 




FIG. 2. The behavior of the adhesive force between two 
rough, "spherical" particles. The three regimes of the force 
vs. volume of the wetting layer are I - Asperity Regime, II 
- Roughness Regime, and III - Spherical Regime. The insets 
show the extent of the wetting region typical of each regime. 



we must set this disjoining pressure equal to the pressure 
inside the contact regions. If the radius of curvature of 
the contact meniscus is r, then 



2Hr 



(24) 



In the asperity regime, t will increase with the meniscus 
radius of curvature. In the roughness regime, however, r 
saturates at a value ~ In, and t will be constant. In this 
regime any added fluid enters the contact region. Finally, 
in the spherical regime, t will again increase. 

Now consider a system with a volume Vi of liquid per 
particle. If we suppose that the spheres are close-packed, 
then each sphere has 12 neighbors, so there is an average 
of 6 contacts per sphere, each with a fluid volume of 
V = Vi/6. The average number of contacts per unit area 
will be (3(j>v /ttR 2 ), where </>y is the volume fraction of 
the particles. Thus for a close-packed lattice, for which 
<f>v = \/2n/6, the adhesive stress will be approximately 



SA 



fA 

V2R 2 



(25) 



and we can now substitute into Eq. ( J13| ) to obtain the 
dependence of the critical angle upon fluid volume. For 
a sandpile of fixed depth D, this is linear in V, 



tan# c 



kTV 



V2l 2 R R 2 pgD 



sec [tan-^fc)] , (26) 



up to a saturation result determined by (see Fig. |2|) 



Now let us consider the experiment of Hornbaker et 
al. They measured the critical angle for a medium com- 
posed of radius 4 x 10~ 2 cm polystyrene spheres with 
small amounts of added oil using the draining crater 
method 0. They found a linear increase in the criti- 
cal angle measured as a function of the volume of added 
oil. They claimed that the failure of their systems was 
at the surface, and concluded that they could account 
for their results by assuming that 99.9% of the fluid was 
outside of the contact zones between the particles. Their 
particles had a surface roughness on the order of lpm. 

We find that the increase in the critical angle is linear 
with the fluid volume in the roughness regime, up to the 
saturation result Eq. (27). We expect the Hamaker con- 
stant for a wetting fluid to be negative, and of the order 
of magnitude of H ~ 10 -20 erg, so no more than a few 
monolayers of fluid should be present on the particle sur- 
faces. Furthermore, in the linear regime, all fluid added 
to the system will enter the particle contacts. Thus we 
disagree with the claim that the overwhelming majority 
of the fluid in this case will be outside the contact zones. 
Finally, we disagree with the interpretation of their ex- 
periment, according to which surface failure is the most 
relevant failure mode- the failure plane in cohesive ma- 
terials should be at depth. 

We are grateful to P. Schiffer for providing us with 
copies of Ref. || before publication and to D. Erta§ for 
many useful discussions. 
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